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Abstract 

The Faddeev equations for the atomic helium-trimer systems are solved 
numerically with high accuracy both for the most sophisticated realistic po- 
tentials available and for simple phenomenological potentials. An efficient 
numerical procedure is described. The large-distance asymptotic behavior, 
crucial for weakly bound three-body systems, is described almost analyt- 
ically for arbitrary potentials. The Efimov effect is especially considered. 
The geometric structures of the bound states are quantitatively investi- 
gated. The accuracy of the schematic models and previous computations 
is comparable, i.e. within 20% for the spatially extended states and within 
40% for the smaller 4 He-trimer ground state. 
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1 Introduction 



Weakly bound two- and three-body systems in low angular momentum states are 
spatially very extended jl], [|. We shall call them halos, if a substantial 
fraction of the wave function is outside the range of the potentials. Each of 
the constituent "particles" may be a composite system, provided it is a tightly 
bound structure with a binding energy large compared to the interparticle binding 
energy. Then the corresponding (intrinsic and halo) degrees of freedom decouple 
and a few-body treatment is appropriate. Prominant examples are the new class 
of nuclear states called nuclear halos || |J and various molecular systems, in 
particular the atomic helium trimer || |9|]. 

Halo states are in general characterized by their large spatial extension. They 
could as well appear as excited states of systems which in the ground state are of 
a different nature and spatially confined. Borromean systems, where no binary 
subsystem is bound, are particularly interesting three-body halo candidates [|l(| . 
They have by definition a relatively low binding energy. Many nuclear examples 
are available [@, [5], |6| and until recently it was unclear whether the helium trimer 
was of Borromean nature ||. The halo concept is general and of interest in 
different subfields of physics H@BISE,@, 0,11,0,1010 E3I- In the simplest case 
of two-body halo systems universal relations between size and binding energy were 
discussed [p], fl3| . Also the general properties of three-body halos are intensively 
discussed [g (I | § 

The halo structure is almost entirely determined by the tail properties of 
the effective interactions between the particles. If the finite-size particles have a 
substantially overlap, the halo degrees of freedom do not decouple and the few- 
body approximations are inaccurate ||, [14]]. The effective two-body interactions 
should then reproduce low-energy scattering properties, i.e. phase shifts or it may 
be sufficient with the correct scattering length and effective range [Q, [14|, [15[] . For 
molecular systems it is possible to compute the two-body effective potentials 
directly from the electromagnetic interaction || [16 . 



The halo structure is technically difficult to compute accurately, since precise 
knowledge of the behavior of the effective radial (three-body) potential at large- 
distance is indispensable and at the same time the small distances, responsible 
for the actual size of the binding energy, must also be properly included. Fortu- 
nately, a method treating the large distances analytically and the short distances 
numerically has recently become available [T3, [T7], [TSL [TjJ. The method is very 



powerful as demonstrated by the succesful investigation of the patological Efi- 
mov effect H, O, p(J. Within the method the Faddeev equations are solved in 



coordinate space in two steps. First the angular part is solved for each hyperra- 
dius, thereby providing the adiabatic basis. This is usually the most difficult part, 
where semianalytic large-distance asymptotic solutions are employed. Afterwards 
the coupled set of radial equations for the expansion coefficients is solved. 

Another somewhat similar method was recently used to study the helium 
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trimer ||. This method exploit the fact that the interaction (in addition to two 
distance-coordinates) only depends on one angle, since the three Euler angles 
describe rotations of the system as a whole. In principle the same adiabatic basis 
is then obtained. However, the computations only include the lowest angular 
basis state which already seems to involve extensive calculations. 

The 4 He-dimer and trimer systems have recently been experimentally estab- 
lished as weakly bound systems although without direct measurement of the 
binding energies pll . However, the size of the 4 He-dimer was found to be exceed- 



ingly large compared to the range of the two-body potential. In the conceptually 
simple experiment transmission was measured through a screen consisting of holes 



with radii comparable to the dimer size p2 |. The binding energy of the dimer 



must consequently be exceedingly small. However, this does not constrain the 



binding energy of the three-body system which, due to the Thomas-effect [[23 
still can be much larger than the binding energy of the two-body system. On the 
other hand, results from solving the Schrodinger equation as well as momentum 
space Faddeev calculations strongly suggest that the trimer binding energy in- 
deed is very small ||, 23]. Both computations conclude that the trimer has two 



bound states of which the excited state has the characteristica of an Efimov state: 
When the attraction is increased it disappears into the continuum away from the 
discrete three-body spectrum. 

Interesting systems can be formed by substituting one or more of the 4 He- 
atoms by 3 He-atoms. The electron structure is identical and the effective two- 
body interaction is then unchanged. Due to the smaller masses the three-body 
states move up in energy and calculations predict that 4 He 3 He 2 and the 3 He- 
trimer both are unbound ||. There is no simple relation between the size and 
the binding energy of these systems and detailed calculations of the three-body 
wave functions are needed. 

The purpose of the present paper is three-fold. First, we shall repete the in- 
vestigation of Esry et al. || and extent it to include sufficiently many adiabatic 
basis states to obtain precise solutions. We shall keep an eye on the possible 
occurrence of Efimov states. Second, we shall define simple two-body potentials 
reproducing the low-energy scattering behavior of the realistic potentials. We 
shall assess the accuracy of such approximation schemes, which reduce the nec- 
essary computational efforts by orders of magnitude. Third, we shall describe an 
improved method able to deal efficiently with large repulsive cores as traditionally 
used in the parametrization of the two-body molecular interaction. 

After the introduction we describe first briefly in section 2 the general method 
of solving the Faddeev equations using hyperspherical coordinates. The angular 
equations and the effective radial potential at large distances are then analyzed in 
sections 3 and 4. The rather delicate numerical procedure and the corresponding 
results are discussed in sections 5 and 6, respectively. Finally section 7 contains 
the summary and the conclusions. 
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2 Theoretical method 



We shall consider an angular momentum zero system of three interacting inert 
"particles" , where the fc'th particle has mass m k and coordinate r^. Their intrinsic 
degrees of freedom are frozen and only the three-body degrees of freedom shall 
be treated here. The two-body interactions between the particles % and j are V^-. 
We shall use the three sets of Jacobi coordinates Xj and and the corresponding 
three sets of hyperspherical coordinates (p, fli) = (p, ai, fl Xi , Q Vi ) defined by P, |lOf 

■rrijTj + m k r k 
w ' w ' rrij + rrik 

Xi = p sin ttj , yi = p cos ttj , (3) 

/l rn^-mfc \ 1/2 /l 77ij(mj +m fc ) \ 1/2 
Pi* = T ' = , 4 

where k} is a cyclic permutation of {1, 2, 3} and /i 2 are the reduced masses 
of the subsystems in units of a normalization mass m. The volume element in 
terms of one of the sets of hyperspherical coordinates is given by p 5 dQdp, where 
dfl = sin 2 a cos 2 adadQ x dQ y . The kinetic energy operator is 

n, & ( -5/2 d 2 5/2 15 A 2 \ 

T = ^\- p ' W p + ' (5) 

1 r) 2 P P 

A 2 = — t^tt^ sin(2a) + -J§- + — |- - 4 , (6) 
sin (2a;) acr sin a cos^ a 

where the angular momentum operators P x and I 2 are related to the x and y 
degrees of freedom. 

The total wave function \l/ of the three-body system is written as a sum of 
three components which in turn for each p are expanded on a complete set 
of generalized angular functions 3>W(/>, 

* = E ^)( Xl , y 4 ) = 4j £ /„(p) E $«(p, ^) = -4 E fi ) . ( 7 ) 

i=l P n i=l P n 

where p _5//2 is the radial phase space factor and n = 1, 2, 3... 

The angular functions, the adiabatic basis, are now for each p chosen as the 
eigenfunctions of the angular part of the Faddeev equations: 



h 2 



\ (A 2 - \ n (p)) *co + Vjk (®w + *w + $(*)) = o , (8) 



2m p 

where k} is a cyclic permutation of {1, 2, 3}. 
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The radial expansion coefficients f n (p) are obtained from a coupled set of 
"radial" differential equations 0, i.e. 



dp- 



d 2 2m(E-V 3 (p)) If 15\ \ 

+ -J (A»(p) + —) -Qnnj fn{p) 
= £ ( 2Pnn>-^- + Qnn) fn'(p) , (9) 



where E is the three-body energy, V^(p) is an anticipated additional three-body 
potential and the functions P and Q are defined as angular integrals: 



P nn ,(p) = J d^ n (p^)— $ n ,(p, fi) 



o2 

Q w (p) = | dn$*(p,fi) — ^fofi) • 



<9p 2 



(10) 



(11) 



Explicit expressions for these coupling terms can be obtained by deriving 
ec f(H) w ith respect to p. By use of the normalization condition ($„|$„') = 8 nn t. 
we obtain 



Pnn i Pnn' 



d( P 2 v) 
dp 



$n' 



for n n' 



Qnn ^ ] PnmPmn j 
m^n 

\ ^ A n A m 
w = 7 _ ^ 

m^n,m^n' n n 



PnmPmn' 



9 2 p 



$n') + 2 ($ 



<9p 



<9(p 2 ?;) 



$n' ) Pnn' 



An — A n ' 



(12) 
(13) 

(14) 

for n n' , 



where v essentially is the sum of the three two-body potentials, i.e. v = Hf=i v i 
with Vi(x) = 2mVjk(x/ pjk) /h 2 . Thus P is antisymmetric. 



3 Large-distance angular solution 

The angular functions $^(p, Qi) are expanded in products of the three-body 
spherical harmonics Yi xmx {Q, Xi ) and Y£ ymy (Q yi ) where the orbital angular mo- 
menta and their projections associated with x and y are (£ x ,m x ) and (£ y ,m y ). 
The total angular momentum is zero and we get 

^\p, «i) = E ^frrv ® W^)] 00 » ( 15 ) 
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where sin(2aj) is a factor related to phase space. We can now insertion of eq.([T5|) 
into eq.(||) and obtain an equation for each of £ projections of each Faddeev 
component: 

/ d 2 £(1+1) £(£ + 1) 2 , . 2/ A , , . 

+ — ji — + 2 — +P Wsinai) -v n (/>) ]<p)^{p,oti)= (16) 



\ daf sin a>i cos 2 y 

-p^psinaA-R^aA , 

where i?^ is the sum of the other two Faddeev components projected onto 
\Yt m (Q Xi ) (g) Yt- m (Q yi )] 00 in the i'th Jacobyan system. I.e. i?^ is a sum of all 
the non-zero £ -contributions of the two other Faddeev components. A detailed 
analysis shows that R^\ai) is proportional to sin^ +1 c^. 

When p is large the potential is zero everywhere except for small a«. Thus 
for £ > the effect of the potential on the left hand side of the equation have to 
been seen through a centrifugal barrier proportional to £(£ + 1) and the effect of 
the potential on the right hand side is quenched by a factor of a\ +l . 

Thus for non-zero £ all terms containing the potential can be put to zero 
for large enough p. Thus the free equation appear. So if any of the Faddeev 
components have a finite contribution with £ > for p — > oo we obtain the free 
solution for that component with the free eigenvalue v 2 n = (K + 2) 2 where K is 
an even integer equal or greater than 2£. Again a detailed analyses shows that 
v\ = (K+2) 2 + 0(p~ 1 ~ 2£ ). This analysis is actually quite similar to the ones given 
below except higher transcendental functions have to be used in the derivation. 

If we now assume that all the Faddeev components only have contributions 
from £ = for p — > oo we end up with three coupled asymptotic equations 

/ d 2 \ 

- p 2 Vi(p sin a>i) + is 2 (p) <f$(p,ai) = p 2 Vi(p sin a { ) 



1 fn-/2-\ir/2-<p k -ai\ , 



tanyjfc 



\ 



fTlk{ m i + m j + Wife) 



miTrij 



The short-range potential Uj(p sin c>A can be assumed to vanish for on larger than 
a constant «q which for large distances must decrease inversely proportional to 
p. The corresponding solution to eq.(|l7|) is then 



<ffi'(p,aii) = sin 



7T 



for ai > af , (19) 



where Aq is an arbitrary constants and the proper boundary condition </>q (p, a 
7r/2) = explicitly is included. 
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The function 0g (p, ay) in eq.(|T7|) is only needed when fj(psinaj) remains 
finite, i.e. for on < ckq . When o$ + aig < (fk the integration in eq. fll7D only 
includes the region where eq.flTTD is valid. Therefore as soon as p corresponds to 
distances outside the short-range potentials we can rewrite eq.([l7|) as 

/ d 2 \ 

p 2 Vi(p sin an) + v 2 (p) <fio \p,ai) = 2p 2 Vi(p sin 



xE^ r " v . v; ■ (20) 



da 2 r -^--v • wy-ru vr,-v ~r -, vr ~--. y p 

■ sin (y((p k - §)) 
j^i sm(2<f k ) 

The solution to eq. (p0|) is given in eq.(^) for ctj > and for < we 
find 

d) a) 2 A^-* sin (^(^fc — |)) 
4>o (p, oii) = (f>$l(p, at) - - sin(va i ) ^ sin(2y? fc ) ' ^ 

where is the solution to the homogeneous equation 

( d 2 \ i 

- p 2 t>;(p sin o^) + u 2 (p) ^(p, a { ) = (22) 
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with the boundary conditions 4>qI{p, cxi — f ) — 0o/!(P' a « = 0) = 0. If fj(psinaj) 
varies within the interval a» < qjq the solution to eq.(|2"2]) could be complicated 
and only obtained numerically. However, if Vi(x) approaches a finite constant for 
small x we find immediately 

41 (p, a,) = 4° sin (a^v 2 {p)-p 2 vm) (23) 

for arbitrary constants -Bo*''- 

Assuming that vct^ ' is small for large p we can expand the right hand side of 
eq.flTT]) to leading order in ctj. This gives 

(Jaf _ p 2 v i (psm.a i ) + ^ 2 (p)^j (j>o\p, ad = 2ajC (i) p 2 ^(p sin a { ) , (24) 

^(i) _ 4>o\p,Vk) + 0Q fc) (p, ^j) , 25 x 

sin(2<£ fc ) sin^j) 

The eigenvalue solutions v 2 to eqs.([[7|), ( ^0|) and (|24|) converge towards the hyper- 
spherical spectrum as p increases. Due to the coupling the asymptotic values are 
now approached over a distance defined by the scattering lengths, which might 
be very much larger than the ranges of the interactions. 

The potentials p 2 fj(p sin a-j) vanish for large p for all a-j except in a narrow 
region around zero. The conditions for the effective range approximation therefore 
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become better and better fulfilled as p increases and any potential with the same 
scattering length and effective range would lead to the same results. Let us then 
in the region of large p use square well potentials Vjk(r) = — So ©(r < Ri), or 
equivalently expressed by the reduced quantities Vi(x) = —s^Q(x < X{ = RiPjk), 
where the parameters are adjusted to reproduce the two-body scattering lengths 
and effective ranges of the initial potential. The corresponding solutions are then 
accurate approximations to our original problem at distances comparable to the 
range of the potentials |L9] . 



Using square well potentials we can now easily define the constants 0$ = 
arcsin(Xj/p) such that the potentials Vj(psinaj) for large p are zero in region II 
where a« > 0$ and finite and constant in region I where ccj < 0$ <C 1. The 
square well solutions are then given in eqs.([T9|), (pT]) and fl23[). We can also solve 
the approximate equation in eq. (f24"D for both region I and II. We obtain again 
the same solutions, where v 2 {p) is neglected compared to p 2 t>,(0) and sin(i/«j) in 
the last term of eq. (|21~l) should be replaced by voti. Matching these solutions at 
the boundary a, = olq provides the quantization condition and determines the 
eigenvalue v. 

Also the decoupled equations in eq . (|i~6D for non-zero I- values can be solved 
for square well potentials. In region I and II the solutions are proportional to the 
spherical Bessel functions and the hyperspherical polynomials, respectively. 

All these square well solutions coincide at large distances with the solutions 
for any general potential with the same scattering length and effective range. 
However, when the two-body potentials are strongly repulsive at small distances it 
is an advantage to use eq. (|2tj|) with the correct potential instead of both the square 



well solution and the additional approximation in eq. (|24|) . A faster convergence 
(for smaller p) towards the large-distance asymptotics is then obtained. 

The method described here is useful when the large-distance behavior is im- 
portant. In extreme cases as for example for Efimov states and halo systems the 
large distances are necessary and otherwise difficult or impossible to compute 
sufficiently accurate. For large distances, where the approximations become in- 
creasingly better, it is in any case much easier to obtain accurate solutions from 
the asymptotic equations than from the original integro-differential equations in 
eq. (|T7|) . Thus we have formulated a convenient numerical procedure. 



4 The effective potential and the Efimov states 

The most important ingredient in the radial equation is the angular eigenvalue 
A = v 2 — 4. We can extract the essential properties by analytical analysis of a 
schematic model. For simplicity we consider three identical bosons. We first find 
the square well potential —SoQ(r < i? ) with the same scattering length a s and 
effective range R e ff as the original two-body potential. The lowest-lying adiabatic 
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radial potentials are then asymptotically close for the two potentials for hyperradii 
p larger than 2Rq. The large-distance behavior is essentially determined by the 
s- waves. The resulting eigenvalue equation, derived from eqs.([HJ), (^T|) and 



is a simple transcendental equation [19 



- cos{a K) 



sm{a Q K) 



sin(z/ 7r/6) sin(j/a ) — v sin ((^o — vr/2)t- 
^ sin(z/ 7r/6) cos(a ^) — ^cos ((a — tt/2)z/) , (26) 



where k = \J2mS p 2 /fi 2 + z/ 2 (p) and a = arcsin(i? /pv / 2)- 

The lowest-lying hyperradial potential in eq.(^) are to the lowest orders in 
1/p obtained by expansion of eq . (|26|) . For Rq <^ p <^ \p,ya s \, where a s is the 
scattering length, we find 
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ce - ^ 



A(p)«i4-4-c B ^(l-^), (27) 

p mv E 

z/|sin(^z/ £ ) / N 

— r — E / 2 f 7 - = 0.486214 , 28 

^ cos( ~ u E ) - cos(f i/ s ) + \v E sm(^v E ) 

where z/g = —1.0125 is the solution to the Efimov equation obtained by expansion 
of eq.(^6"D in the limit where a s /p = ±oo, i.e. 

8sin(z/ s vr/6) = u e V3cos{u e tc/2) . (29) 

In the other limit where \pijd s \ <C p we find instead for the lowest eigenvalue[] 

_ _ A8a sf j, ij for ^ > Q A _ p^ for ^ < Q (30) 

vrp ' p^af 

The effective radial potential is in this way obtained analytically for the square 
well potential for various distance intervals when the scattering length is large 
compared to the interaction range. However, as indicated above these results 
are more general, since they can be expressed in terms of scattering length and 
effective range, which are universal quantities characterizing low-energy two-body 
scattering properties of short-range potentials. For distances smaller than Rq the 
details of the original two-body potentials are essential and the square well results 
are not applicable. This part of the potential is decisive for the actual (small) 
binding energy of the three-body system. 

Between a few times Rq and a few times below 15pij\a s \ eq.(|2"7]) is valid. This 
region could be called the Efimov region, since the Efimov states appear here 
for p < \pijd s \ where the radial potential in eq.(0) behaves as p~ 2 . The lowest 



1 We assume a negative scattering length for a two-body system with one bound state 
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angular and radial wave functions $ and f n are then obtained from eqs.fll^D, 
and (PI) as 



3 



sm(v E (aci - 7r/2)) 



$««E sm(2^) ' /n K ^ Kv ' {PV 2m \ E \/ h2 ) > ( 31 ) 

where K UE is the modified Bessel function with an imaginary index. The radial 
function could easily be obtained from eq.(^) in the limits where either the en- 
ergy or the centrifugal barrier term is dominating. The same result is found by 
expansion of eq.(|31~D for small and large values of p 



f n oc ^sin ( \v E \ ln( P - ) J for p 2 < h 2 /(2m\E\) . (32) 



eff 



f n oc exp(-p^2m\E\/h 2 ) for p l > fr/(2m\E\) , (33) 

where the zero point for the first oscillation in p is assumed to be pijR e ff. The 
exponential fall off at large distance occurs for all bound states for distances 
larger than h 2 / (2m\E\), since the effective radial potential falls off faster than 
>> 

The number of Efimov states is determined by the number of possible oscilla- 
tions of the sine function within the interval p < |15//^a s |. Each new oscillation 
corresponds to another further excited state and the zero point for the fc'th os- 
cillation is then given by 



p k w pijReffexp yj^J ■ (34) 
The number of Efimov states N is estimated from eq. (|30D by pn ~ |15//ya s |, i.e. 

iV^lnf^. (35) 



R 



eff I 

The fc'th state has a radial extension of about p^, i.e. exponentially increasing 
with k, and correspondingly an exponentially decreasing energy 

h 2 ( 2kn\ , . 

2l7lR eff V We\J 

Clearly infinitely large a s then gives infinitely many bound Efimov states. They 
are extremly weakly bound and extremely spatially extended. 

Efimov states may also be present in non-identical bosonic three-body sys- 
tems. The condition is that at least two of the binary subsystems have bound 
s-states at zero energy or equivalently infinitely large scattering lengths. The 
third binary subsystem must be without bound states. A finite number of states 
with similar properties appear when the conditions are approximately fulfilled. 
The appearance is determined by the large-distance behavior of the zero angular 
momentum states. 
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5 Numerical procedure 



It should be straightforward to solve the equations of motion numerically. How- 
ever, two major difficulties must be tackled. The first is related to the large 
short-range repulsive core where the wave function must be extremely small. 
The potential has an attractive pocket at short distances where the wave func- 
tion also is small. These inner regions are where the potential is relatively large, 
are decisive for the energy of the bound states and the corresponding wave func- 
tion must therefore be calculated with high accuracy. The second problem is that 
the weakly bound states are located far outside the range of the potential and 
the appropriate long-range part of both the potential and the wave function are 
therefore also decisive. A marginal numerical inaccuracy could be crucial for the 
number of loosely bound states. 



5.1 Method 

Our first overall computing strategy is to use the Faddeev decomposition of the 
Schrodinger equation. The proper boundary conditions and the more subtle 
correlations are then much easier to account for in actual computations. The 
second computing strategy is to expand on a suitable angular basis for short 
distances and solve the large- distance asymptotic angular equations directly at 
large distances. In an overlapping region these solutions are then combined. 
Finally the coupled set of radial equations is solved. 

For short distances each Faddeev component is then expanded on the hyper- 
spherical basis [|l^, [1^, [|. This amounts for each p to an expansion in eq. (|T5|) 



of the functions </>W in Jacobi polynomials. The remaining expansion coefficients 
are functions of p. The total angular basis then consists of a carefully selected set 
of hyperspherical harmonics with the total angular momentum zero. The basis 
is characterized by the Faddeev components i = 1,2,3, the quantum numbers 
I = 0, 1, 2, 3... and the hyperspherical quantum number K — 0, 2, 4, 6.... 

Only basis functions with the correct symmetry properties are included, i.e. 
symmetric or antisymmetric for bosons or fermions. When p 2 v i is small for small 
p the angular eigenvalue equation eq. (H) only includes the kinetic energy operator 
with the corresponding hyperspherical harmonics as solutions. In this case a few 
of the lowest t and K quantum numbers are sufficient in the basis. 

For a repulsive core, where p 2 Vi diverges (or simply is very large), larger l- 
values are needed in the basis to produce a vanishing wave function at small 
distances. The angular matrix elements of the potentials become large for small 
p even for high £, where the contribution usually vanishes due to the centrifugal 
barrier. A faster convergence in I is achieved by subtracting a constant Vq from 
each of the two-body potentials and subsequently correcting the energy eigenvalue 
by 3Vq. This is simply a shift of the energy scale by 3V . If the infinite basis is 
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used the choice of constant makes no difference. A minimal basis size is obtained 
when V approximately equals the lowest angular eigenvalue. 

Increasing p normally decreases the required maximum £- value since the cen- 
trifugal barrier then increases. On the other hand larger hyperspherical quantum 
numbers (or polynomial orders) are needed, because for short-range interactions 
the potentials in a-space become increasingly confined to a narrow region of phase 
space around ai — 0. Accurate wave functions in this decreasing region around 
ai = require high-order Jacoby polynomials (high K) to reproduce the rapid 
change caused by the corresponding behavior of the potential. 

In the limit of very large p the hyperspherical spectrum is also approached for 
short-range potentials. Therefore in this extreme limit only small £ and small K 
are needed to obtain the angular eigenvalue spectrum. However, in the sometimes 
very extended transition region of intermediate p-values the eigenvalues are re- 
quired with high precision. Thus, increasing p demands an increasing basis for an 
accurate description. It is crucial that this basis only needs to include high-order 
polynomials, but still only small angular momenta. The necessary basis can then 
be substantially reduced from the full hyperspherical basis with all K and £ due 
to the use of the Faddeev equations instead of the Schrodinger equation. 

We shall treat large p by considering s- waves, £ = 0, without the further 
expansion on the Jacobi-polynomia, i.e. we shall solve the large-distance equa- 
tions (f2y) directly. Moving towards smaller p we must then find an overlap 
region where the large- and small-distance solutions coincide. To optimize com- 
putational speed and accuracy the basis must therefore be carefully selected to 
provide sufficient accuracy in the overlap region. 

When the asymptotic states correspond to £ ^ the convergence is expected 
to be much faster and essentially reached already at small distance. We extrap- 
olate the A-values by using the form K(K + 4) + cp~ 2 where c is a constant 
obtained by matching and K is the asymptotic hyperspherical quantum number 
of this eigenvalue. 

The direct way of calculating the coupling terms P nn i and Q nn i is to perform 
a numerical derivation, but this numerical double derivative is inaccurate in re- 
gions where the wave functions change rapidly, i.e. for repulsive cores at small 
distances. Instead we used the explicit expressions in eq.(|T2|) for small distances. 
For large distances we used numerical derivation for the P nn i and Q nn and a 
smooth extrapolation of the small distance functions for the non-diagonal Q nn '- 



5.2 Parameters 



In the actual cases of the He-trimer systems we choose the mass unit m in eq. (|5|) 
in the definition of the Jacobian coordinates 10, to be one nuclear mass 
unit m = 1822.887 a.u. (atomic mass unit). The masses of the He-atoms are 
then m( 4 He) = 4.002603250m and m( 3 He) = 3.016026m. 
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The two-body interaction is the LM2M2 potential from p6[. The three-body 



interaction V3 is very small |3]J and assumed to be zero in the present computa- 
tion. For the 4 He-trimer the effective range of the potential is R e ff = 13.843 a.u. 
(1 atomic length units = 0.529177 A) and the binding energy is estimated to be 
of the order of lmK [g [22| (1 Kelvin = 0.0861735 meV = 3.16679 -lO" 6 atomic 
energy units). This extremely low binding energy is equivalent to a two-body 
scattering length of a s = —189.054 a.u., which is much larger than the effective 
range. The size of the 4 He-dimer is measured to be (|ri — r 2 |) = 120±20 a.u. [22]. 
For the 4 He — 3 He system the different masses instead give a s = 33.261 a.u. and 
Reff = 18.564 a.u. and for 3 He - 3 He we get a s = 13.520 a.u. and R eff = 25.717 
a.u. 

For the 4 He-trimer we used a basis set of 150, 40 Jacobi polynomials for each 
Faddeev component for £ = 0,2, respectively and 30 for £ = 4, 6, 8. This amounts 
to 840 basis states in total, but after projecting on the states with correct 
symmetry the basis size was reduced to 235. For the 3 He 4 He2 and 4 He 3 He2- 
molecules we used a basis set of 150,60,50,45,40,35,30,25 Jacobi polynomials 
for each Faddeev component for i = 0,1,2,3,4,5,6,7, respectively and 20 for 
£ = 8, 9, 10, 11, 12, 13, 14. The components related to two identical 4 He or 3 He- 
atoms only needs even £ due to the symmetry requirement. For the asymmetric 
trimers this amounts to a total basis of 1500 which reduced to 638 after symme- 
try restoration. In all cases this corresponds to an overlap region between small 
and large distances around p = 100 a.u. The converged results from the basis 
diagonalization and the asymptotic solutions are then used below and above 100 
a.u., respectively. 

The convergence and the accuracy of the numerical results were tested thor- 
oughly by varying the number of Jacobi polynomials in each component as well 
as the number of partial waves. The results are also independent of relatively 
small variations of the position of the matching point between the small and large 
distance solutions. Thus we expect reliable numerical results. 



6 Numerical results 

The numerical procedure and the input parameters are now specified. The results 
are obtained by expansion on the adiabatic basis and we shall therefore first 
compute the corresponding angular eigenvalues. Then we shall show energies and 
wave functions for both the LM2M2 and the simple schematic potentials. Finally 
we discuss the various intrinsic geometric structures with special emphasis on the 
possible appearence of Efimov states. 
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6.1 Angular eigenvalues 

The decisive quantities are the adiabatic potentials where the main components 
are the angular eigenvalues. The lowest of these are computed and shown in fig. [I] 
s functions of p for the 4 He-trimer. Both the general small-distance numerical 
solutions and the asymptotic large-distance s-wave solutions are shown. They 
diverge for small p due to the strongly repulsive core in the interaction between 
two He-atoms. The lowest of these eigenvalues has an attractive pocket at small 
distance and diverges parabolically towards — oo for large p. The pocket is re- 
sponsible for the three-body bound states and the divergence reflects the bound 
state in the two-body subsystem. Due to the very small two-body binding en- 
ergy the lowest level remains almost constant up to about 1000 a.u. where the 
divergence sets in. The overlap region above which the asymptotic solutions can 
be used is seen to be around 100 a.u. 

The higher lying eigenvalues asymptotically approach the hyperspherical spec- 
trum defined by K(K + 4) where K is a non-negative integer. In the present case 
K must furthermore be an even number because the total orbital angular mo- 
mentum is zero. The requirement of totally symmetric wave functions removes 
a number of otherwise possible degenerate states. This eliminates completely all 
K = 2 states whereas the asymptotic states corresponding to K = 0, 4, 6, 8, 10 all 
are non-degenerate and with the angular quantum number i = 0. For K = 12 the 
s-wave solution is accompanied by a higher angular momentum solution, which 
has converged to the asymptotic value 12(12+4)=192 already at about p = 100 
a.u. This level and that of K = 10 avoided crossing at about p = 80 a.u. The 
even higher lying levels are also degenerate corresponding to more than one to- 
tally symmetric state. 

The angular eigenvalue spectrum changes when the asymmetric 3 He 4 He2- 
molecule is considered, see fig. |2]. The interactions are the same as for the 
symmetric case but the masses are different. The number of levels is roughly 
increased by a factor of two due to the less restrictive symmetry requirement. 
The divergence for p = 0, the parabolic divergence at large p due to the two- 
body bound state and the convergence to the hyperspherical spectrum at large 
distances are still the dominating features. Compared to the 4 He-system we now 
find a number of additional levels. They correspond to nearly antisymmetric so- 
lutions in the exchange of 3 He and 4 He whereas the "original" levels also present 
for the 4 He-trimer must be of symmetric character. 

Considering the 4 He 3 He2-molecule obtained by substituting another 4 He by a 
3 He-atom again decreases all the binding energies. We show the angular eigen- 
value spectrum in fig. 0, where the usual features are seen. However, now the 
lowest eigenvalue converges to zero for large p, since none of the two-body sub- 
systems form bound states. The asymptotic spectrum is now reached much faster 
than for the two previous cases. This rate of convergence is roughly proportional 
to the sum of the three scattering lengths or more precisely — Z)i=i,3 a i*Vjfc 
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Figure 1: The lowest angular eigenvalues A n for the 4 He-trimer as functions 
of p. The logarithmic scale on the p-axis uses the atomic length unit. The 
long-dashed curves are the angular eigenvalues obtained by diagonalization for 
the relatively small distances, the dotted curves are the large-distance s-wave 
solution and the short dashed curve is the extrapolation of the 7'th eigenvalue 
(non-zero t) obtained by using the form 12(12 + 4) + cp~ 2 . The hyperspherical 
quantum numbers K are shown to indicate the asymptotic spectrum K(K + 4) 
for p — > oo. The 7'th and 8'th eigenvalues both approach the same asymptotic 



value 12(12 + 4) = 192. The lowest eigenvalue diverges as 
the bound dimer state. 
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Figure 2: The lowest angular eigenvalues \ n for the 3 He 4 He2-trimer as functions 
of p. The logarithmic scale on the p-axis uses the atomic length unit. The 
long-dashed curves are the angular eigenvalues obtained by diagonalization for 
the relatively small distances, the dotted curves are the large-distance s-wave 
solution and the short dashed curve is the extrapolation of the 9'th eigenvalue 
(non-zero £) obtained by using the form 8(8 + 4) + cp~ 2 . The hyperspherical 
quantum numbers K are shown to indicate the asymptotic spectrum K{K + 4) 
for p — > oo. The lowest eigenvalue diverges as — p 2 corresponding to the bound 
dimer state. 
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Figure 3: The lowest angular eigenvalues A n for the 4 He 3 He2-trimer as functions 
of p. The logarithmic scale on the p-axis uses the atomic length unit. The 
long-dashed curves are the angular eigenvalues obtained by diagonalization for 
the relatively small distances, the dotted curves are the large-distance s-wave 
solution and the short dashed curve is the extrapolation of the 9'th eigenvalue 
(non-zero £) obtained by using the form 8(8 + 4) + cp~ 2 . The hyperspherical 
quantum numbers K are shown to indicate the asymptotic spectrum K(K + 4) 
for p — > oo. The lowest eigenvalue remains finite at large p. 



is the scattering length of the two-body system with particles j and k), i.e. -2042 
a.u., -459 a.u. and 195 a.u., respectively. The different rates are clearly seen in 
figs. [I], Q and |3] especially when the range of the plotted A- values are considered. 
The levels all approach their asymptotic values from below, since the sum of scat- 
tering lengths in this case is positive, compare to eq. (|30|) for identical particles 
for s- waves. 



6.2 Energies and wave functions 

For halo states the decisive properties of the potentials are related to the low- 
energy scattering behavior or equivalently to the large distance behavior. This 
is formally expressed in the effective range theory where disparate potentials 
give the same low-energy results provided they have the same scattering length 
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Table 1: Strength and range arameters for gaussian, V(r) = — So exp(— r 2 /b 2 ), 
exponential V(r) = —Soexp(—r/b) and square well —SqQ(t < b) poten- 
tials for the different He two-body systems. They all have the same scat- 
tering length a s and effective range R e ff as the LM2M2 interaction [po], i.e. 
a s = -189.054,33.261,13.520 a.u. and R eff = 18.564, 13.843, 25.717 a.u. for 
4 He- 4 He, 4 He- 3 He, 3 He- 3 He, respectively. The strengths are in Kelvin (1 K = 
0.0861735 meV = 3.16679 -10~ 6 a.u.) and the ranges are in atomic units (1 a.u 
= 0.524177 A). 



Potential 


4 He 


- 4 He 


4 He 


- 3 He 


3 He- 


3 He 




So(K) 


b(a.u.) 


So(K) 


b(a.u.) 


So(K) 


b(a.u.) 


Gaussian 


1.227 


10.03 


0.8925 


10.55 


0.16186 


11.27 


Exponential 


3.909 


4.117 


3.284 


4.088 


2.499 


4.233 


Square well 


0.5578 


14.28 


0.3691 


15.61 


0.2399 


17.07 



and effective range. Since smooth and slowly varying potentials allow easy and 
accurate computations we adjusted the parameters of a square well — 5*oO(r < 
Rq), a gaussian V(r) = — S exp(— r 2 /6 2 ) and an exponential potential V(r) = 
— S exp(— r/b) to reproduce a s and R e jf of the LM2M2-potential. The resulting 
parameters are given in table [I]. 

The lowest and decisive eigenvalue for these potentials is plotted in fig. £| as 
function of p. The eigenvalues for all the three simple potentials coincide for 
p > 40 a.u. whereas that of LM2M2 lies slightly below and slowly approaching 
the others as p increases. The deviations are enhanced due to the focus on 
one A- value and the restriction to distances less than 200 a.u. The differences 
can only be interpreted as effects of terms of higher order than in the effective 
range expansion. They are not completely negligible in this region because of 
the large repulsive core of the LM2M2-potential. At smaller distances the simple 
potentials produce eigenvalues approaching zero for p = 0. This behavior is then 
compensated by less attraction in the pocket region. 

The binding energies are obtained by solving the coupled set of radial equa- 
tions. The lower and upper limits can be found by using only the lowest adiabatic 
potential with and without the diagonal term Q u , see f3(J. These limits were 
computed in ||. We repete this computation using our method and we further- 
more increase the number of adiabatic potentials until the results have converged. 
For the 4 He-trimer we also find one excited state in agreement with the discussion 
in connection with eq . (|35|) . One additional oscillation would extend beyond 10 
times the scattering length and a third excited state is therefore not possible even 
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Figure 4: The lowest angular eigenvalue as a function of p for the 4 He-trimer for 
the realistic LM2M2 [^6j potential and the three schematic potentials defined in 
table |I| The square well results are obtained from eq . fl26|) . 
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Table 2: The energies in mK of the ground state and the excited state of the 
4 He-trimer for one realistic and two schematic potentials. The rows in the table 
display the numerical results for different numbers of adiabatic channels in the 
basis. We also give the results of Esry et al. || calculated for the LM2M2 
potential |p6| . The row labeled Qn = displays the results obtained with the 



lowest adiabatic potential without the diagonal term Qn- 



Potential 


Esry et.al §] 


LM2M2 |§ 


Gaussian 


Exponential 


On = o 


-293.7 


-3.518 


-293.5 


-3.517 


-152.6 


-2.576 


-176.0 


-2.823 


l 


-106.1 


-2.118 


-105.9 


-2.121 


-150.2 


-2.468 


-173.9 


-2.714 


2 






-119.5 


-2.224 


-150.6 


-2.482 


-174.1 


-2.729 


3 






-122.0 


-2.245 


-150.7 


-2.485 


-174.2 


-2.731 


4 






-123.9 


-2.259 


-150.7 


-2.485 


-174.2 


-2.731 


5 






-124.3 


-2.263 


-150.7 


-2.485 


-174.2 


.2.732 


6 






-124.7 


-2.265 


-150.7 


-2.486 


-174.2 


-2.732 


7 

8 






-125.1 
-125.2 


-2.267 
-2.269 











without the Qn-term. 

We compare in table ^| binding energies for the LM2M2 and the simple inter- 
actions for both the ground state and the excited state of the 4 He-trimer. We 
confirm the limits computed in |§. The lowest 7 to 8 adiabatic potentials are 
needed to obtain the accurate energies deviating from the upper limit by about 
20% and 7% for the ground and excited states, respectively. This relatively slow 
convergence for the LM2M2 potential is due to the large repulsive core. Less adi- 
abatic potentials are needed for convergence for the simple two-body interactions, 
where the lower and upper limits in any case are very close and a larger basis 
therefore also for this reason is less needed. The binding energies for the gaussian 
potential are about 20% and 10% higher than for the LM2M2 interaction for the 
ground state and the excited state, respectively. The corresponding numbers for 
the exponential potential are 40% and 20%. 

The large- distance tails of gaussian and exponential two-body potentials are 
very different and the variations established here therefore measure the uncer- 
tainties inherent in such estimates. These results strongly indicate that such 
weakly bound states can be studied by simple but carefully chosen potentials 
provided an accuracy of less than about 40% and 20% are required for states lo- 
cated respectively at the edge and far outside the short-range potential. Further 
accuracy can be obtained with an appropriate choice of radial shape of the simple 
two-body potential. The computations simplify due to the lack of repulsion at 
small distances and the smaller basis required. 
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Table 3: The energies in mK of the ground state of the 3 He 4 He2-trimer for one 
realistic and two schematic potentials. The rows in the table display the numerical 
results for different numbers of adiabatic channels in the basis. We also give the 
results of Esry et al. || calculated for the LM2M2 potential p6| . The row labeled 
Qu = displays the results obtained with the lowest adiabatic potential without 
the diagonal term Qu- 



Potential 


Esry et.al [|] 


LM2M2 |p|l 


Gaussian 


Exponential 


Qu = 


-86.52 


-86.47 


-20.44 


-26.06 


l 


-10.22 


-9.682 


-18.41 


-24.27 


2 




-10.12 


-18.67 


-24.52 


3 




-11.90 


-18.74 


-24.57 


4 




-12.74 


-18.79 


-24.60 


5 




-12.86 


-18.81 


-24.61 


6 




-13.23 


-18.81 


-24.62 


7 




-13.34 


-18.81 


-24.61 


8 




-13.55 


-18.82 


-24.61 


9 




-13.64 


-18.82 


-24.61 


10 




-13.66 


-18.82 


-24.61 



Substituting one of the 4 He-atoms in the 4 He-trimer by a 3 He-atom decreases 
the binding energies and only the ground state remains bound. We compare in 
table [3] the results for different potentials. Unlike the symmetric case we now 
obtain slightly less binding energy than in ||. The difference is less than 0.1% 
compared to the strengths of the potentials but about 5% compared to the binding 
energy itself. Our convergence was thoroughly tested again, but we were not able 
to find the origin of the discrepancy. As for the symmetric system an accurate 
binding energy requires 5 to 6 adiabatic potentials for the LM2M2 interaction 
and again about 40% is gained compared to the upper limit obtained by using 
only the lowest adiabatic potential. The simple potentials are now off by factors 
1.4 and 1.8 for gaussian and exponential radial shapes, respectively. 

The radial wave function has a distribution of components corresponding to 
the different adiabatic potentials, see eq.(0). The lowest is by far carrying the 
largest probability as seen in tables f| and |5| for the two He-trimer systems. The 
exceedingly small occupation probability for the higher lying angular eigenvalues 
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Table 4: Occupation probabilities as functions of adiabatic channels for the 
ground state (left) and the excited state (right) of the 4 He-trimer for the three 
different potentials in table |]. 



Potential 



LM2M2 H| 



Gaussian 



Exponential 



1 
2 
3 
4 
5 
6 
7 



0.9978 
1.72-KT 3 
4.0-10~ 4 

8.5- KT 5 
2.4-KT 5 

4.6- KT 6 
1.0-10~ 5 
3.9-10" 6 



0.9989 

1.02- 10 -3 
2.740" 5 
5.540" 6 
2.640" 6 

8.3- 10- 7 
1.2-10- 7 
3.240" 7 



0.9997 
2.61-10" 4 

2.8- 10" 5 
1.5-10~ 6 

5.9- 10~ 7 
1.4-10~ 7 



0.9991 
9.16-10" 4 
1.940" 5 
4.3-10~ 6 
2.0-10~ 6 
5.840- 7 



0.9998 
1.69-10- 3 

1.6- 10~ 4 
8.9-10~ 7 

3.7- 10~ 6 
9.4-10- 8 



0.9990 
9.840- 4 
1.840" 4 
3.740" 5 
8.440" 7 
2.640" 7 



seems off hand to contradict the relatively large contribution to the binding en- 
ergies. The rather instructive explanation for weakly bound systems is that the 
variation in binding energy might appear to be large when compared to the bind- 
ing energy itself, but insignificant compared to the strength of the interactions 
responsible for the binding, see fig. |] where the potential depths are around 2 
K compared to the binding energies of mK. Thus the lowest adiabatic potentials 
give a good approximation for the wave functions although the corresponding 
binding energies only are obtained with less relative accuracy. 

The dominating components of the radial wave functions for the ground state 
and excited state of the 4 He-trimer are shown in fig. |5] along with the different 
effective radial potentials. The repulsive core pushes the ground state towards 
larger distances, whereas the simple potentials allow finite probability also at 
small distances. The excited state is already at a large distance and therefore less 
affected by the changing potentials. The wave functions related to the different 
potentials are quite similar. The length of the first oscillation is about the 350 
a.u. as predicted from eq. (|34]) with R e ff = 18 a.u. The next would extend 20 
times as far, i.e. about 7000 a.u. 



6.3 Geometric structure 

The wave functions in fig. [| are clearly not revealing directly the geometric struc- 
ture of these three-body systems. This is not only because the angular part of the 
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Table 5: Occupation probabilities as functions of adiabatic channels for the 
ground state of the 3 He 4 He2-trimer for the three different potentials in table 
i 



Potential 


LM2M2 


Gaussian 


Exponential 


1 


0.9959 


0.9974 


0.9979 


2 


3.43-10~ 3 


2.50-10" 3 


2.03-10" 3 


3 


2.8-10" 4 


1.0-10" 4 


7.5-10~ 5 


4 


2.4-10~ 4 


2.6-10~ 5 


1.2-10~ 5 


5 


6.9-10~ 5 


1.6-10" 5 


9.9-10~ 6 


6 


4.7-10~ 5 


5.4-10~ 7 


3.0-10~ 7 


7 


1.9-10- 6 


8.9-10- 7 


5.7-10- 7 


8 


1.9-10" 5 


1.2-10" 6 


6.9-10~ 7 


9 


2.4-10" 6 






10 


6.1-10" 7 
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Figure 5: The lowest adiabatic effective radial potentials (thick curves) V e ff(p) = 
fi ((Ai + f)p- 2 ~ Qn) and the wave functions corresponding to the two bound 
states (thin curves) for the 4 He-trimer for one relalistic |26| and two schematic 
two-body interactions defined in table |I| The scale on the p-axis is logarithmic. 
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Figure 6: The contour diagrams of the particle density distributions in the in- 
trinsic coordinate system for the ground state (left) and the excitede state of the 
4 He-trimer (middle) and the bound state of 3 He 4 He2 (right). The units on the 
axes are atomic lengths units. They are different for the three states. There is a 
factor of 1.5 between each contour curve. The scales on the plots are different. 

wave function is missing, but also because the plotted total wave function fully 
includes all the appropriate symmetries, i.e. symmetric in exchanges of identical 
bosons and rotational symmetry averaging over all directions. To get a visual 
picture of the structure we want to consider the plane through the three particles 
at points r^. The center of mass is at the origin. Let us take the y-axis as the 
principle axis with the lowest moment of inertia. We furthermore assume that 
one particle is in the first and two particles are in the third and fourth quadrants. 
This leaves us with the intrinsic geometry of the trimer. 

The resulting contour plots are shown in fig. || Note that the distance scales 
on the individual figures are quite different. The triangular shape for the 4 He- 
trimer ground state (left) is quite apparent. The excited state in the middle show 
that two of the particles prefer to be close while the third particle is further away. 
The Efimov conditions are almost fulfilled, i.e. very large scattering lengths for at 
least two of the binary subsystems. As we shall discuss later this figure illustrates 
the structure of an Efimov state. Any higher-lying Efimov state could simply be 
found by a scaling of the axes, see eq. (|31|) . 

The contour diagram at the right side of fig. |6| also show the preference for 
two relatively close-lying particles. However, only one scattering length is now 
large and the Efimov conditions are not fulfilled. It is instructive to imagine that 
the 3 He- 4 He scattering length is increased and the subsequent appearence of the 
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Figure 7: The contour diagram of the 3 He density (left) and half of the 4 He 
density distribution (right) in the intrinsic coordinate system for the bound state 
of 3 He 4 He2 (right). The units on the axes are atomic lengths. There is a factor 
of 1.5 between each contour. Both the scale and the contours are the same as on 
the total density distribution in figure |[ 



Efimov states. The lowest eigenvalue in fig. |2| would still bend over and diverge 
due to the bound state. The second eigenvalue would approach a negative con- 
stant above the lowest eigenvalue. The corresponding adiabatic potential would 
isolated provide the Efimov states. However, these states would appear in the 
continuum as excited states above the ground state shown in the figure. Thus the 
ground state of the asymmetric system cannot be an approximate Efimov state. 
The superficial resemblance of the middle and right figures therefore illustrates 
how deceiving the density distribution can be. Additional information is needed 
to come up with the correct interpretation. 

For the asymmetric system we can also compute the individual particle den- 
sity as shown in fig. [7|. The left figure can be obtained from the right figure if 
a small component with an elongated distribution is added. The 4 He-trimer ap- 
pears in a completely symmetric state while the asymmetric system in addition 
also has a relatively small asymmetric component of the 3 He- 4 He relative wave 
function. The large symmetric component resembles the 4 He-trimer and the small 
asymmetric configuration, which becomes dominating at larger distances where 
the structure of the 4 He-dimer is approached. 

The triangular structure at the left side of fig. ^| is apart from the scale similar 
to the three-body cluster approximation of the ground state of 12 C, whereas the 
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two other trimer states are very different from the almost linear configuration of 
the astrophysically interesting first excited + -state of 12 C, which therefore is far 



from being an Efimov state, see [^J . 

The geometry of the He-trimers can be studied quantitatively by computation 
of various expectation values related to the distances between the particles. We 
first consider the sizes shown in table |6| where two measures are used, i.e. the 
root mean square radius and the (first order) distance from the center of mass. 
The dimer has an average distance between the particles of about 100 a.u., which 



is at the lower end of the experimentally allowed range of 120 ± 20 a.u. ||22|| . The 



three-body ground state is much smaller than the 4 He-dimer || 22]. However, 



the size of about 100 a.u. of the excited state is around 50% larger than for the 
dimer both when measured by the first and the second moment. This radius is 
larger than the range of the effective two-body interaction but smaller than the 
size (189 a.u.) of the scattering length. The spatial extension of the probability 
reaches beyond \a s \ as seen in fig. [5], but the wave function is still well within 
the range of the effective radial three-body potential (p~ 2 behavior) extending to 
about |15/Xy-a s |, see eq.(^0|). These proportions are characteristic for an Efimov 
state. 

The ratio between the three-body radii of excited and ground state is 9.5 which 
only is two times smaller than the asymptotic prediction from eq . (|34D . In com- 
parison the ratio of the binding energies is 52 which is about 7 times smaller than 
predicted in eq. (|36|) for two Efimov states. These observations simply confirm 
that the ground state does not have the characteristics of an Efimov state. An 
additional (second) excited state (which is impossible for this scattering length) 
would have obeyed these scaling laws for radii and energies for Efimov states. 

The asymmetric system 3 He 4 He2 has only one bound state with both binding 
energy and size in between the corresponding values for the ground state and the 
excited state of the symmetric system, see tables 0, |3| and ||. For Efimov states 
the product of the energy and the square of the radius should be a constant, see 
eqs. flH^D and fl3T5|) . This product is four times larger for the excited symmetric 
state than for the asymmetric state. The radius should then be about twice as 
large if the state should resemble an Efimov state. However, this seems to be 
unlikely, since the radius is of the same order as the two-body effective ranges. In 
any case the interpretation as an approximate Efimov state was already excluded 
in the discussion above. Thus a pronounced halo structure is therefore a more 
appropriate description. 

The distance between the two 4 He-particles computed as (|ri — r2|) is 28 a.u. 
and the distance between 3 He and one specific of the 4 He-particles is 38 a.u. The 
distances between both the 3 He- 4 He and the two 4 He particles are then two to 
three times smaller in this asymmetric bound state than in the bound 4 He dimer. 
As necessary both distances are in between the minimum and maximum of 20 and 
46 a.u. Again this demonstrates that the bound dimer-state is not a dominating 
configuration in the three-body systems. 
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Table 6: Expectation values of various operators related to the geometrical struc- 
ture of the bound states of the He-dimer and the He-trimers. The results for the 
different potentials in tables and |5| are given for each quantity: LM2M2 (top), 
gaussian (middle) and exponential (bottom). The coordinate of particle i is r i; 
the center of mass coordinate is R and n = 2, 3 is the number of particles. All 
the lengths are in atomic units. 





4 He 2 


4 He 3 


4 He^ 


3 He 4 He 2 




66.9 


11.8 


115 


26 


V / (E i (r i -R)^)/n 


66.9 


10.8 


104 


22 




66.9 


9.9 


95 


20 




49.0 


10.2 


96 


20 




49.0 


9.4 


86 


17.5 




49.0 


8.6 


77 


15.6 




98.1 


11.8 


70 


20 


{vaiatfj \rij\) 


98.0 


10.6 


67 


19.1 




98.0 


9.6 


64 


17.0 




98.1 


23 


219 


46 




98.0 


21 


195 


40 




98.0 


19.3 


178 


35 




1 


0.51 


0.34 


0.46 


/ min^j |ry| \ 
Xmax^j |r<j| / 


1 


0.51 


0.37 


0.49 




1 


0.51 


0.38 


0.49 
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Another interesting quantity in table |6] is the expectation value of the ratio 
between the smallest and the largest interparticle distance. First the minimum 
and maximum average distances themselves differ by a factor of two for the 
ground states and a factor of three for the excited Efimov-like state. The ratio 
is nearly 0.5 for both ground states and 0.34 for the excited state. If the angular 
part (including the a coordinate) of the wave function was constant this number 
would be 0.52. In the case of an Efimov state this number was computed to 
be 0.38. Thus the excited state does also in this respect resemble an Efimov 
state. That the ratio is less than the value for an Efimov state indicates non- 
negligible contributions corresponding to two particles in a bound state and the 
third particle further away. This is consistent with the observation that the 
excited stated is larger than the dimer. 

In all cases the simple potentials produce sizes deviating about 10% and 20% 
for gaussian and exponential potentials from those of the LM2M2 interaction. 
Thus if only this accuracy is required the corresponding simpler computations 
are sufficient. 



7 Summary and conclusions 

Spatially extended three-body halo states are apparently present in some of the 
atomic helium-trimer systems. The 4 He-dimer is experimentally established as 
exceedingly large and as a consequence it must also be very loosely bound. The 
large distances are then essential in any accurate description. The possibility of 
finding Efimov states is adding interest in obtaining a quantitative understanding, 
but unfortunately at the same time emphasizing the practical difficulties. The 
necessary large- distance correlations must be treated very accurately without 
loosing the precision at smaller distances where the binding energy originates. 
This three-body problem is of quantum mechanical nature since the probability 
almost entirely is found in the classically forbidden region. 

We incorporate the correlations by solving the coordinate space Faddeev equa- 
tions. We expand the wave function in a complete set of angular eigenfunctions 
obtained by solving the angular part of the Faddeev equations for a given hy- 
perradius. This adiabatic basis is related to a set of effective radial potentials. 
The corresponding radial wave functions are solutions to the coupled set of radial 
differential equations. We first formulate this prescription in general and then 
the crucial large distances are specifically considered. Only s-waves couple at 
large distances and this fact is exploited to derive simple equations with almost 
analytic solutions. A suitable numerical procedure follows from this formulation. 
For exploratory calculations we furthermore suggest a simplifying approximation 
where the "correct" potential is substituted by almost schematic potentials with 
relatively little loss of accuracy for the spatially extended states. 
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The lowest effective radial potential is decisive and we describe a series of its 
large-distance analytic properties for a symmetric system. When the two-body 
scattering length is sufficiently large the Efimov states appear as solutions. We 
repete the occurrence conditions and a few of the analytic properties of these 
exotic states. 

The strongly repulsive core is necessary in a correct description of the two- 
body potential. The numerical procedure is somewhat delicate and a sizable basis 
is necessary. Each of the Faddeev components is expanded on the hyperspherical 
basis. We need relative angular momentum quantum numbers up to about 10, 
Jacobi polynomials for s-waves up to order 150 and significantly fewer for the 
higher partial waves. This numerical procedure produces the solutions at small 
distances up to about 100 a.u. and from then on we use the solutions obtained 
directly from the decoupled large-distance set of equations. 

The schematic two-body potentials are parametrized to reproduce the scatter- 
ing length and effective range of the original complicated potential. We use radial 
shapes of gaussian, exponential and square well form. The solutions are much 
easier to compute due to the lack of hard core repulsion. Still these simple at- 
tractive potentials reproduce the results of the "correct" calculations within 40% 
and significantly better for the spatially extended states. More suited shapes 
could certainly reduce this uncertainty, which in fact is remarkable for several 
reasons. First this demonstrates the importance of the large-distance behavior of 
the effective radial potential. Second the error in binding energy by using such 
simple potentials is roughly the same as the results obtained by using only the 
lowest adiabatic potential, which carries more than 99% of the probability in a 
fully converged calculation. Thus simple potentials could be exploited to gain 
even semi quantitative insight. 

We have essentially confirmed the calculations of [§] using the lowest adia- 
batic potential. In our full calculations we find additional binding of 7% to 40% 
measured relative to the binding energy, which is several orders of magnitude 
smaller than the strength of the attraction. Indeed these states are very weakly 
bound. We investigated the geometry of the bound He-trimer states. The ground 
state of 4 He 3 is an equal sided triangle. The excited state of 4 He 3 and the ground 
state of the asymmetric 3 He 4 He2 system are both also triangular but with sides 
differing by a factor of two to three. 

The excited state of 4 He3 has all the characteristics of an Efimov state whereas 
the somewhat similar density distribution of the asymmetric state rather should 
be characterized as a pronounced halo state. More Efimov states would only 
appear in the symmetric system when the scattering length is increased by about 
one order of magnitude. Smaller changes would be needed to promote the halo 
state to the giant halo state called an Efimov state. 

In conclusion, we have described a method to compute spatially extended 
three-body structures. The essential large-distance behavior is specially treated. 
We apply the method to the atomic helium trimers. The binding energies are 
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accurately determined and the geometric structures are found to be triangular 
configurations. The ground state of 4 He3 has a triangular shape and the excited 
state is an Efimov state whereas the 3 He 4 He2 bound state is a halo state. 
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